arXiv: 1504.04363vl [math.CO] 16 Apr 2015 


TRANSFINITE FORD-FULKERSON ON A FINITE NETWORK 


SPENCER BACKMAN AND TONY HUYNH 


Abstract. It is well-known that the Eord-Fulkerson algorithm for finding a maximum 
flow in a network need not terminate if we allow the arc capacities to take irrational 
values. Every non-terminating example converges to a limit flow, but this limit flow 
need not be a maximum flow. Hence, one may pass to the limit and begin the algorithm 
again. In this way, we may view the Ford-Fulkerson algorithm as a transfinite algorithm. 

We analyze the transfinite running-time of the Ford-Fulkerson algorithm using ordinal 
numbers, and prove that the worst case running-time is For the lower bound, 

we show that we can model the Euclidean algorithm via Ford-Fulkerson on an auxiliary 
network. By running this example on a pair of incommensurable numbers, we obtain 
a new robust non-terminating example. We then describe how to glue k copies of our 
Euclidean example in parallel to obtain running-time cj*. An upper bound of is 
established via induction on \E\. We conclude by illustrating a close connection to 
transfinite chip-firing as previously investigated by the first author [2]. 


1. Introduction 

The Ford-Fulkerson algorithm m is a classic algorithm for computing the maximum 
flow in a network. At each step, the algorithm hnds an augmenting path P from the 
source vertex s to the sink vertex t and then pushes as much flow as possible along P. 
This procedure is then iterated, until no such augmenting path exists. It is well-known 
that in certain networks with irrational capacities, the Ford-Fulkerson algorithm does 
not necessarily terminate if the augmenting paths are not chosen carefully. The smallest 
non-terminating example is due to Zwick |24| . 

Dinits m and Edmonds and Karp m independently showed that if one always chooses 
an augmenting path of minimum length, then the Ford-Fulkerson algorithm will necessarily 
terminate. We warn the reader that whenever we refer to the Ford-Fulkerson algorithm, 
we mean the original version where augmenting paths can be chosen arbitrarily. 

It is fairly easy to show that if the Ford-Fulkerson algorithm does not terminate, then 
it will converge to a (not necessarily maximum) flow /. Thus, after uj steps, we may begin 
the algorithm anew, starting with the limit flow /. By iterating this procedure, we can 
view the Ford-Fulkerson algorithm as a transhnite algorithm and ask what its worst case 
running-time is in terms of ordinal numbers. Note that the notion of using ordinals as a 
complexity measure dates back at least to the work of Turing |23j . 

The following theorem, which is the main result of this paper, determines this worst 
case ordinal running-time up to a constant factor in the exponent. 

Theorem 1. The worst ease running-time of the Ford-Fulkerson algorithm is 

This theorem is established via the following two lemmas. 

Lemma 2. For every network N, every run of the Ford-Fulkerson algorithm terminates 
after at most steps. 

The proof of Lemma is by induction. Although we are working with ordinal numbers 
where transhnite induction might seem like a natural tool, the argument proceeds by hnite 
induction on the exponent \E{N)\. 
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Lemma 3. For every i € N, there exists a network N{tj on I ares and a run of the 
Ford-Fulkerson algorithm on N{i) with run-time at least cjLasJ. 

To prove Lemma[^ we first construct a new non-terminating example of Ford-Fulkerson. 
The main idea is to demonstrate that the Euclidean algorithm can be modelled by applying 
Ford-Fulkerson to a particular network. Thus, a run of the Euclidean algorithm on two 
incommensurable numbers gives a non-terminating example of Ford-Fulkerson. To obtain 
the general result, we demonstrate how to glue several copies of the Euclidean example in 
parallel and run them lexicographically. 

In |2], the hrst author investigates a certain transhnite chip-hring process on metric 
graphs. In the last section we describe a close connection between the results presented 
in this article and those appearing in [2]. 

2. Ordinal Running-time 

We now describe how we intend to measure running-time via ordinal numbers. For 
an introduction to ordinals and network flow theory, we refer the reader to |18| and [9], 
respectively. Roughly speaking, we use an extended notion of a Turing machine that can 
complete an inhnite number of steps of computation, and continue computing afterwards. 
This matches the notion of infinite time Turing machines by Hamkins and Lewis |16| . 
However, as a tradeoff, we allow a saboteur to choose the augmenting paths at every step 
of the algorithm. The ordinal running-time is then the worst running-time over all possible 
sets of choices of the saboteur. We give the precise details below. 

All networks considered will always be hnite with hnite arc capacities. Each step of 
the Ford-Fulkerson algorithm will be indexed by an ordinal a and the corresponding flow 
after step a will be denoted /„ (we allow /q be any valid flow). A run of the algorithm 
is obtained by choosing an augmenting path Pp+i for fp (if it exists) for each successor 
ordinal (3 + 1, and setting //j+i to be the flow obtained from by pushing as much flow 
as possible along Pp+i- If no augmenting path exists at step /3 -|- 1, we dehne fp+i to be 
fj 3 . If a is a limit ordinal, we dehne /„ to be a certain limit how. Some care must be 
taken to ensure that this limit how is well-dehned and this is the content of Lemma [d] 

The run-time of a particular run is the least ordinal a such that fa = /a+i- For a hxed 
network N, the (worst case) running-time of the Ford-Fulkerson algorithm on N is the 
maximum of the run-times over all runs of the algorithm on N . The main result of this 
section is that transhnite Ford-Fulkerson is a well-dehned procedure. 

Lemma 4. For every ordinal a, every run of the Ford-Fulkerson algorithm assigns a 
well-defined flow fa after step a. 

Proof. Let N = (V, E) be a network with a source s £ R, a sink t € V, and non-negative 
hnite capacities c{u, v) for each arc (u, v) £ E. We proceed by transhnite induction. We 
let /o be any valid initial how. Now, let a be an ordinal, and assume that fp is dehned 
for all ordinals (3 < a. 

First suppose a is a successor ordinal, say a = (5 If there is no augmenting path for 
fjS, then we set fa ■= fg- Otherwise, if the saboteur chooses the augmenting path 
for fg, then we dehne fa to be the how obtained from fg by pushing as much as possible 
along P/ 3 + 1 . 

If a is a limit ordinal, we proceed as follows. For each arc e and ordinal fi < a, we 
dehne a { — 1 , 0 , l}-valued variable y{e, fi) as follows. If /3 := 7 -|- 1 and the saboteur chose 
Py+i as the augmenting path for /+, then we set y(e, fi) to be 1 if e is a forward arc of P/ 3 , 
-1 if e is a backward arc of P/ 3 , and 0 otherwise. We initialize y{e,0) := 0. We also let xq 
be the value of the initial how /o and set X /3 to be the amount of how pushed by Pg. If /3 
is a limit ordinal or there is no augmenting path at step /3, we set both y(e, fi) and xg to 
be 0 . 
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Consider Yli3<a^^- Observe that at most countably many terms are non-zero, since 
this sum is bounded (by the capacity of a minimum cut). Moreover, this series converges 
absolutely, since each term is non-negative. Therefore, this sum is independent of the 
order of summation and is hence well-defined. We define a flow /„ by setting 

fa{e) = /o(e) + ^y(e,/3)x^ 

l3<a 

for each e € E. Observe that X]/3<a absolutely convergent series, since 

converges. Therefore, X]/ 3 <q 2/(^)/^)®/3 i® ^1®° ^ well-defined sum. Evidently, 
0 < fa{u,v) < c{u,v) for all {u,v) £ E since 0 < ffs{u,v) < c{u,v) for all jS < a. It 
remains to verify that /„ satisfies conservation of flow. Let u £V \ {s,t}. Observe that 

V fa{v,u)- V fa{u,v)= V lim/^(r;,R)- V lim/^(u,r;) 

p—¥a p^a 

{v,u)£E {u,v)£E {v,u)£E (u,v)gE 

= Yj //3(^W)) 

{v,u)^E {u,v)^E 

= 0 

where the last equality follows from conservation of flow for /^. □ 

3. The Upperbound 

Given a flow /, an arc e is a zero-arc if /(e) = 0, is saturated if /(e) is equal to the 
capacity of e, and is extreme if it is saturated or a zero-arc. 

Lemma 5. For every k £ N, every network N, and every run of the Ford-Fulkerson 
algorithm on N, either the algorithm has already terminated after steps or there are at 
least k -£ 1 extreme arcs in f^k 

Proof. We proceed by (finite) induction on k. Since some arc must be extreme after 
pushing as much flow along an augmenting path, the lemma clearly holds for k = 0 (note 
= 1). We inductively assume that the claim holds for k. Now, for each j £ N, let 
Aj be the set of extreme arcs for the flow fj^^k. Since ui^ steps have passed between 
ju!^ and (j -|- l)a;^, by induction we may assume that \Aj\ > k -£ 1 for all j. Next note 
that each e £ A can only switch between being a zero-arc and a saturated arc a finite 
number of times because the value of the flow is always bounded. This implies that for 
all sufficiently large j, there is an arc aj ^ Aj such that Oj is extreme in fj^kj^i- Let A 
be such that Aj = A infinitely often and a be such that Oj = a infinitely often. Since 
lim^_,.^^fc+i fjs = limj^t^; fj^k, it follows that each edge e £ A\s extreme for f^^k+i. On the 
other hand, since /^fe+i = linij^i^ fj^k_^i, we conclude that a is also extreme for f^k+i. 
Thus, fi^k+i has at least k -\-2 extreme arcs, as required. □ 

Using Lemma we now prove our upperbound, restated for convenience. 

Lemma For every network N, every run of the Ford-Fulkerson algorithm terminates 
after at most steps. 

Proof. If the algorithm has not terminated after steps, then Lemma [^implies that 

for all j £ N, every arc of N is extreme in /j^|B(iv)i-i. Let c be the smallest non-zero 
capacity over all arcs of N. Since all arcs are extreme in fji^\E(N)\-i, the flow increases 
by at least c at step -£ 1 (for all j). This is a contradiction since the value of 

fuj\EiN)\ is finite. □ 
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4. The Lowerbound 

Let a,b € M_|_. We begin by constructing a network Na^b and a run of the Ford-Fulkerson 
algorithm on Na^b which simulates the Euclidean algorithm on a and b. 

s s 




Figure 1. A network Na^ and two augmenting paths which allow for 
the subtraction of b from a. 


Lemma 6. Let Nafi be the network depicted in Figure^ U ^ ^ Qj then there is a run of 
the Ford-Fulkerson algorithm on Na^ with run-time at least io. 

Proof. The labelled arcs of Figure (denoted Ca and Cb) have capacities a and b with 
0 < b < a. The other arcs have large capacity (which we will specify later). Given a flow 
/ and an arc e, the residual capacity of e, denoted c'{e), is defined to be c(e) — /(e). Let 
Pi and P 2 be the augmenting paths shown in red. Now, starting from the zero-flow, if we 
push flow along Pi and P 2 , then c'{ea) becomes a — b and c'{eb) is still b. We continue this 
process until ai := a — nib < b. At this point, note that the value of the currrent flow is 
2nib and that the roles of ai and b are reversed. 

Let P{ and P 2 be the reflections of Pi and P 2 through the vertical line from s to t. 
By next pushing flow along the augmenting paths P{ and P^, we can convert the residual 
capacity of Cb to b — ai. We continue this process until 61 ;= b — miai < ai. Note that the 
horizontal arcs are backward arcs of P{ and P^. However, since 2miai <2b< 2nib, there 
is enough flow along the horizontal arcs to perform these 2mi augmentations. The roles of 
ai and bi have been reversed again and we continue inductively. Therefore, this example 
“computes” the greatest common divisor of a and b. If | ^ Q, it will have run-time at 
least cj, as required. However, it remains to check that the total flow is still bounded after 
u steps (so that we may specify the capacities). 

Claim 1. The total flow after uj steps of the Euclidean run on Na^b is o-t most 4(a -|- b). 

Proof. Let a = ao > b = bo > ai > bi > ... be the intermediate outputs of the Euclidean 
algorithm. Observe that for all n > 0, On+i < and bn-ei < Therefore, the total 
flow after ui steps is at most 2 Yl’^=oi^n + + ^)) required. □ 

Thus, for all arcs e ^ {cq, Cb}, we can take c(e) = 4(a -|- b). □ 

This example is quite robust in that a small perturbation of the arc capacities will almost 
certainly produce another non-terminating example. We will now “glue” k copies of Na^b 
in parallel to obtain a run of the Ford-Fulkerson algorithm with run-time at least This 
proves Lemma A key property of our gluing construction is that we can “recharge” one 
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copy of the Euclidean example using only two steps of the Euclidean algorithm of another 
copy. 

Lemma For every ^ G N, there exists a network N{i) on i ares and a run of the 
Ford-Fulkerson algorithm on N{i) with run-time at least . 

Proof. For each /c G N, we define a network and show that there is a run of the 
Ford-Fulkerson algorithm on with run-time at least uj^. For simplicity, we only 
illustrate the case k = 2] the general case is similar. Figure depicts the corresponding 
network for k = 2. The labelled arcs have capacities a, b, a, and b respectively with b < a 
incommensurate. We denote these arcs as £i,..., from left to right. The remaining arcs 
have large capacities, which we will specify later. We denote the three middle vertical arcs 
as 61,62 and 63 from left to right. 

Note that there are “left” and “right” copies of Na^ sitting in depicted in green 

and blue, respectively, in Figure]^ 


s 




Figure 2. A network (depicted twice) that admits a run of Ford- 
Fulkerson that takes iterations. Left and right copies of Na^ in are 
shown in green and blue. 

Let a = ao > b = bo > ai > bi >... be the intermediate outputs of the Euclidean 
algorithm for (a, b). We begin by running the Euclidean algorithm on the left copy of Nafi. 
This takes oj steps. We then run two steps of the Euclidean algorithm on the right copy of 
Aq b, obtaining residual capacities ai > bi for £3 and £ 4 , respectively. We next perform a 
recharging step (described below), and then re-run the (full) Euclidean algorithm on the 
left copy of Na^b- We then run another two steps of the Euclidean algorithm on the right 
copy of Nafi, obtaining residual capacities 02 > 62 and proceed iteratively. 
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The Recharging Step. At the beginning of the nth recharging step the residual 
capacities of and £4 are 0 , 0 , a„ and bn, respectively. We hrst nse the two 

angmenting paths in red to recharge the residual capacity of ii withont drastically 
changing the rest of the network. Observe that none of the middle arcs 61,62 and 63 
are in the left or right copy of Nafi- The arcs 62 and 63 will always have zero-flow except 
dnring this recharging step. By pnshing flow along the two red angmenting paths, c'(ei) 
becomes and c'(62), c'(e^) and 6^(64) are nnchanged. Note that both 62 and 63 have 
zero-flow after the recharging step, as claimed. Similarly, nsing two angmenting paths 
similar to the ones in red, we can recharge c'(62) to bn withont changing c'(ei), c'(63), and 
c '(64). 

One easily checks that there is always snfhcient flow along backward arcs to rnn the 
Enclidean algorithm on the right copy of Na^b and to perform the recharging steps. 
Fnrthermore, since a and b are incommensnrate, Un and bn are also incommensurate for 
all n. Therefore, this rnn of Ford-Fnlkerson on reqnires at least uP steps. By Claim 

the total flow after steps is at most 

00 00 

4 (a -|- 6 ) + 4 'y ^ (ctn + bn) + 2 ^ + bn) < 14 (a -|- b), 

n=0 n=l 

where the hrst term corresponds to the Enclidean rnn on the right copy of Nafi, the 
second term corresponds to all the Enclidean rnns on the seqnence of left networks 
^aofio, ■ ■ ■, and the third term corresponds to the total how prodnced by the 

recharging steps. 

We can therefore take c(e) = 14 (a -|- b) for all e ^ {ii,... ,£4}. The proof is complete 
as it is straightforward to check that contains at most 25 /c edges for all k. □ 

Note that onr example for the lowerbonnd starts with the zero-how, while onr proof for 
the npperbonnd is valid starting with any initial how. 

5. Chip-firing and a Restricted Duality 

The hrst anthor | 2 ] proved the non-termination, and investigated the transhnite 
behavior, of a certain greedy chip-hring algorithm on metric graphs. Many of the resnlts 
obtained and the argnments employed in that paper are similar to those which appear in 
this article. In this section we attempt to shed some light on the similarity between 
these two works by showing how in a restricted setting, a variant of the chip-hring 
algorithm previonsly investigated is planar dnal to the Ford-Fnlkerson algorithm. We 
wish to emphasize that the dnality presented in this section is not strong enongh to imply 
that the resnlts of either paper follow from the other. We begin with a review of chip-hring. 
We have attempted to make onr introdnction elementary and historically motivated, bnt 
nearly all proofs are omitted. There are many references provided for the interested reader. 

Chip-hring is a simple game played on the vertices of a graph, which has been 
independently discovered in several different commnnities la El Ea EI]. A configuration 
of chips on the vertices of a graph is a fnnction from the vertices to the integers. We 
imagine each integer as describing some nnmber of poker chips sitting at the corresponding 
vertex. The vertices with a negative nnmber of chips are said to be in debt. Given a chip 
conhgnration D, a vertex v fires by sending a chip to each of its neighbors and losing its 
degree nnmber of chips in the process so that the total nnmber of chips is conserved. This 
setnp gives a combinatorial langnage for nnderstanding the integer translates of the lattice 
generated by the colnmns of the Laplacian matrix. There has been a recent explosion of 
interest in chip-hring as the natnral langnage for developing a theory of divisors on graphs 
urn and abstract tropical cnrves caEniiizj analogons to the classic theory for algebraic 
cnrves. An abstract tropical cnrve is essentially a metric graph; an edge weighted graph 
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Figure 3. A sequence of two augmenting paths that recharges the far left 
vertical edge is shown in red. 

where each edge is isometric to an interval of length equal to the associated weight. A 
chip configuration on a metric graph is a function from the points of the metric graph to 
the integers with finite support. In the discrete case, if a set of vertices U fires, a chip is 
sent across each edge in (U, U^). For a metric graph, we define the basic chip-firing moves 
by taking an edge cut of length e and sending a chip across each segment in the cut. 

In keeping with the language of (tropical) algebraic geometry, we may refer to a chip 
configuration as a divisor, and we say that two divisors are linearly equivalent if we can 
get from one to the other by a sequence of chip-firing moves. The main combinatorial tool 
for studying chip-firing on metric graphs is the q-reduced divisor, whose discrete analogue 
is also known as a G-parking function | 22 ) and is dual to the recurrent configurations in 
the sandpile model 0 0 . We say that the number of chips in a divisor D is the degree of 
D and write deg(Z)) for this quantity. 

Definition 1 . A divisor D is g-reduced if D{p) > 0 for allp q and for all A C V{G)\q, 
firing the set A causes some vertex to he sent into debt. 

The fundamental theorem about g-reduced divisors is the following. 

Theorem 7 . | 5 ] |T 5 ] [T 7 ] [ 2 n| Each divisor D is linearly equivalent to a unique q-reduced 
divisor. Moreover, D is linearly equivalent to a nonnegative divisor if and only if the 
equivalent q-reduced divisor is nonnegative. 

In the discrete setting, Dhar nni provided an efficient method for testing whether a 
divisor is g-reduced. Imagine that each vertex has D{v) firefighters present. A fire is 
started at the root q which spreads through the graph, but is prevented from burning 
through the vertices by the firefighters present. When the fire approaches a vertex v from 
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more than D{v) directions, the firefighters are overpowered and the fire passes through 
the vertex. Dhar observed that D is reduced if and only if the fire consumes the entire 
graph. If D is not reduced, then the set of vertices which are not burnt can be fired (not 
to be confused with burned) simultaneously without sending any vertex into debt, thus 
bringing the divisor closer to being reduced. It is easy to check that if a divisor D is 
not reduced and Ui and U 2 are two sets of vertices which can be fired without sending 
any vertex into debt, then we can also fire Ui U U 2 without sending any vertex into debt. 
Thus, the set of “fireable” vertices forms a join semilattice, and Dhar’s algorithm finds the 
unique maximum element in this semilattice. 

Given a divisor on a metric graph which is not g-reduced, one may repeatedly perform 
maximal firings towards q. If this process terminates, we arrive at the unique equivalent 
g-reduced divisor. Unlike the discrete case, it is not clear a priori whether this process 
will terminate. Luo |19| introduced a metric version of Dhar’s algorithm and showed that 
this gives a finite method for computing the associated g'-reduced divisor. The question 
of whether the greedy reduction method also terminates in finite time was left open. The 
first author P] demonstrated that the greedy reduction algorithm need not terminate, but 
as is the case with the Ford-Fulkerson algorithm, the greedy reduction algorithm always 
has a well-defined limit. Thus, the reduction algorithm was interpreted as a transfinite 
algorithm and its running time was analyzed using ordinal numbers. The main result of 
[2] is the following. 

Theorem 8. |2] The worst case running time for the greedy reduction algorithm of a 
divisor on a metric graph is _ 

Given an orientation O of a discrete graph, we can associate a chip configuration Do 
by taking the indegree minus one at each vertex. Our connection between Ford-Fulkerson 
and the greedy reduction algorithm is derived from the following well-known relationship 
between g-reduced divisors and acyclic graph orientations. 

Theorem 9. A divisor D with D{q) = —1 is a maximal q-reduced divisor if and only if 
D = Dq where O is an acyclic orientation with a unique source at q. 

If we take a directed cut {U, U‘^) in O and reverse all of the edges, we obtain a new 
orientation O' such that the associated chip configuration Dq' is obtained from Dq by 
firing all of the vertices in U. Every acyclic orientation is equivalent via source reversals 
(or more generally directed cut reversals) to a unique orientation with a unique source 
at g. This object can be obtained by arbitrarily flipping sources other than g (or more 
generally directed cuts orientated towards g). 

It is easy to check that an acyclic orientation O has a unique source g if and only if every 
other vertex is reachable from g by a directed path. This suggests the following efficient 
method to obtain the unique g-connected orientation equivalent to O by cut reversals. We 
first perform directed search from g. If each vertex is reachable by some directed path, we 
are done. Otherwise, let U be the set of vertices reachable from g. The edge cut (U, U'') 
is oriented towards g, and we can flip this directed cut, bringing us closer to the desired 
orientation. As is the case with chip configurations, given an orientation O and a root g, 
we can associate the collection of sets of vertices C/, such that ([/, U'') is a directed cut in 
O with g G U'', and this collection forms a join semilattice. The cut which we obtain by 
search corresponds to the unique maximum element in this join semilattice. In this sense, 
the search-based algorithm just described may naturally be viewed as a variant of Dhar’s 
burning algorithm for graph orientations. For the remainder of this section we will call 
this algorithm the pseudo Dhar’s algorithm, which is not to be confused with the oriented 
Dhar’s algorithm appearing in [3] (see Section 7 of that paper for a different connection 
between chip-firing and network flows). 
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So as to not mislead the reader, we also highlight some differences between Dhar’s 
algorithm and the pseudo Dhar’s algorithm. The following distinctions hold for both 
the discrete and metric versions of these two algorithms. The first difference is Dhar’s 
algorithm proceeds by firing sets of vertices not including g, whereas the cut reversals 
described in the oriented version correspond to borrowings (the inverse of firings) by sets 
of vertices not including q. The second difference is that Dhar’s algorithm for reducing 
divisors presumes that the divisor in question is nonnegative away from g, but for acyclic 
orientations, the associated divisor is not effective away from q. In fact, the oriented variant 
of Dhar’s algorithm applied to the acyclic orientation produces an acyclic orientation with 
an associated divisor which is nonnegative away from q precisely when q is the unique 
sink, and this occurs when the divisor is g-reduced. 

A metric graph (orientation) can be considered as a limit of (orientations of) discrete 
graphs under repeated subdivisions of the edges, where the ratios of the induced path 
lengths converge to the desired ratios of lengths. On the other hand, networks with real 
capacities may be viewed as a limit of networks having all edge capacities 1 where we 
continue adding parallel directed edges so that the ratios of the number of parallel edges 
converge to the desired ratios of the capacities. For plane graphs, these limits are planar 
dual as replacing an edge with k parallel edges is dual to replacing the dual edge with a 
path of length k. See Figure]^ 



Figure 4 . The duality between planar parallel directed edges and planar 
induced paths of directed edges 


We note that while the termination of Luo’s metric version of Dhar’s algorithm is not 
obvious, the termination of the metric version of the pseudo Dhar’s algorithm is. Every 
time we are stuck with some cut oriented towards q and we push this cut back as far 
as we can, we reach some new vertex. Hence, this process will terminate in at most \V\ 
steps. It is now natural to ask (analogous to the greedy reduction algorithm), if we reverse 
maximal cuts directed towards q arbitrarily to obtain the unique equivalent g-connected 
acyclic orientation, whether this process will terminate in finite time. We will call this 
method the greedy eut reversals algorithm. 

A s-t planar network is a planar network such that s and t belong to the same face. 
We will call an acyclic s-t planar network such that every edge belongs to a directed path 
from s to t a simple s-t planar network. Note that the networks appearing in our lower 
bound construction are of this form. Given a planar network N and an orientation of the 
plane, we define N to be the planar dual network. 

Berge investigated s-t planar networks and showed that a variant of Ford-Fulkerson 
exists which necessarily terminates in finite time and has better worst case running¬ 
time than the Edmonds-Karp variant of Ford-Fulkerson. His method was elegant; when 
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performing the Ford-Fulkerson algorithm, always choose the uppermost augmenting path. 
We are now ready to describe the main observation of this section. 

Main Observation 1. Given a simple s-t planar network N, add an are e = (t, s) of very 
large capacity which passes above N. The Ford-Fulkerson algorithm on N is dual to the 
greedy cut reversal algorithm on the metric graph orientation N U {e}. Moreover, Berge’s 
algorithm for N is planar dual to the pseudo Dhar’s algorithm for N U {n}. 

We note that there is some ambiguity in this duality because the capacity of e has 
not been specihed. However, since v was assumed to have very large capacity, the dual 
directed edge will never affect a run of the pseudo Dhar’s algorithm. Hence, the dual 
of N is essentially well-dehned. By adding {e} to N, we obtain an auxiliary directed 
graph which is strongly connected. Now, augmenting flows in our original network can be 
extended via {e} to reversals of weighted directed cycles containing the unique clockwise 
face q containing s and t. Thus, the general Ford-Fulkerson algorithm may be interpreted 
as greedily performing directed weighted cycle reversals until all of the cycles enclosing 
q are oriented counterclockwise. The planar dual N U {e} is an acyclic metric graph 
orientation and the weighted directed cycle reversals in N (which necessarily contain e) 
are planar dual to directed cut reversals in N U {e} toward q. One may check that there 
are no clockwise directed cycles in our original digraph if and only if the planar dual of 
the auxiliary network is g-connected, where as an abuse of notation we have identihed the 
face q and the corresponding dual vertex. Berge’s uppermost augmenting path extends 
via {e} to a directed cycle which is dual to the cut given by the pseudo Dhar’s algorithm. 
By taking the symmetric difference of this orientation with our original auxiliary directed 
graph and deleting {e}, we obtain a maximum flow in our original network. 

As an application of this duality, by taking the planar dual of the networks and 
augmenting paths described in Section we obtain acyclic orientations such that the 
greedy cut reversals algorithm takes at least steps to terminate. See Figure 



Figure 5 . The ayclic metric orientation dual to the Euclidean network. 
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